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The vortex contribution to the dc field (H) dependent microwave surface impedance Z s = 
R s + iX s of YiBa2Cu307_a thin films was measured using suspended patterned resonators. Z S (H) 
is shown to be a direct measure of the flux density B(H) enabling a very precise test of models of flux 
penetration. Three regimes of field-dependent behavior were observed: (1) Initial flux penetration 
occurs on very low field scales Hi(4:.2K) ~ 100 Oe, (2) At moderate fields the flux penetration into 
the virgin state is in excellent agreement with calculations based upon the field-induced Bean critical 
state for thin film geometry, parametrized by a field scale H S (4:.2K) ~ J c d ~ 0.5 T, (3) for very 
high fields H ^> H s , the flux density is uniform and the measurements enable direct determination 
of vortex parameters such as pinning force constants a p and vortex viscosity r\. However hysteresis 
loops are in disagreement with the thin film Bean model, and instead are governed by the low field 
scale Hi, rather than by H s . Geometric barriers are insufficient to account for the observed results. 



PACS: 74.60.Ge, 74.60.-w, 74.60.Jg, 74.76.Bz 



> ! I. INTRODUCTION 

O 

The dynamics of vortices is important over a wide spectrum of time-scales (or frequencies) . On very long time scales 
Q\ , the problem is relevant to dissipation at or near the critical current. On the other end of the spectrum, at short time 
scales or high frequencies, an understanding of vortex dynamics is important both for fundamental understanding 
and for practical applications of microwave devices. 

Yet another need is to understand flux profiles and field penetration into thin films. Recently there have been 
several calculations based upon critical state models which have yielded field and current distributions in thin films 
for magnetic fields applied perpendicular to the films, and have shown that there are qualitative differences from bulk 
geometries. While some of these results have been tested in actual experiments Jl|, not all aspects of these models are 
fully understood. Most experimental studies of flux penetration and the critical state have involved measurements of 
magnetic properties such as the dc magnetization, and there are few investigations of the dynamic properties. 
c/2 ■ This paper reports on a novel approach, utilizing microwave surface impedance measurements in patterned mi- 
crostrip resonators, to understanding vortex dynamics and flux penetration. Specially designed microstrip resonators 
? *~j ■ patterned from thin films are used as high Q structures to probe vortex response at microwave frequencies in the 
presence of dc magnetic fields. The high sensitivity of these structures enables measurement of the vortex response 
in fields from < O.lOe to several T. The complex microwave surface impedance Z s = R s + iX s was measured as 
functions of magnetic field H, temperature T and frequency u (in the range 1 to 20 GHz). Z s is a measure of the 
total flux in the sample, and hence can lead to precise tests of flux penetration in thin films. In addition the results 
also yield information on vortex parameters such as pinning force constants and viscosity. 

In this paper we focus on results for the changes in surface impedance induced in YiBa2Cu307_5 thin films as the 
externally applied magnetic field is slowly increased from zero, as well as the hysteresis loops which ensue when the 
field is subsequently reduced. The response clearly shows evidence of two relevant field scales. As the field is increased 
an initial linear rise gives way at Hi to a super-linear behavior followed above H s by a linear rise which persists to 
the highest fields available. While virgin response for H > Hi and the magnitude of H s appear to be consistent with 
recent analytical expressions for thin strips in perpendicular fields, the hysteretic response is not so easily understood. 
The critical state model predicts a counter-clockwise hysteresis loop governed by a field scale of order H s , while the 
experiments observe a clockwise hysteresis loop governed by the much smaller field scale H{. Geometrical barriers 
are insufficient to explain the observed hysteresis loops. The results suggest that pinning at the film edges may be 
weaker than in the bulk. 
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II. EXPERIMENTAL DETAILS 



The experiments were designed specifically to enable high sensitivity measurements in dc magnetic fields from 
< lOe to several T. The essential element is the suspended line resonator housed in a Cu package. Microwaves are 
coupled in inductively, by means of loops at the ends of coaxial lines. Similar resonators were used in studies of 
high field || and non-linear || response, which have already been published. We note several key features of the 
experiment: 

• The absence of superconducting ground planes (such as would be present in a traditional stripline or microstrip 
resonator). This allows us to observe effects even at very low applied magnetic fields. 

• The presence of several modes enables measurements at several frequencies in the 1 — 20 GHz range. 

• The ability to carry out both linear and non-linear response experiments independently, and also combined 
rf+dc field experiments where the rf and dc field scales are comparable. 

• Variable coupling at all temperatures enabling weak coupling over a broad range of resonator Q 's. 

Experiments were initially carried out in fields of 1 mT < /io-ff < 4 T in a Janis SuperVaritemp Helium flow 
cryostat. The magnetic fields were generated using a NbTi coil (capable of fields up to 6 T) which is permanently 
mounted in the cryostat. Current was supplied to the solenoid by a Cryomagnetics IPS-50 Magnet Power Supply, 
under computer control via RS-232 using a Cryomagnetics CIM-16. Magnet calibration runs using a LakeShore 450 
Gaussmeter and a 5 ft. cryogenic axial Hall probe showed that residual fields of the order of 10 — 100 Oe/A may be 
present at the center of the superconducting NbTi magnet once it had been energized. For experiments at true zero 
applied static field, and for greater sensitivity at low fields the probe was adapted to fit in a similar flow cryostat 
built by Cryo Industries of America. In this series of low-field experiments, the magnetic fields were supplied by a 
liquid N2 cooled Cu solenoid capable of fields up to 400 Oe when supplied using an HP 6024A power supply. The 
CIM-16 interface module was also used here to control the supply current remotely via RS-232. Current (and hence 
field) reversal was accomplished by manually interchanging the magnet leads at zero current. The field-to-current 
ratio was determined using a LakeShore 450 Gaussmeter and an axial probe to be 106 ±1 Oe/A 

A number of additional features were incorporated into the probe in order to increase the system's sensitivity at 
low magnetic fields. A Minco thin film heater with a nominal resistance of 50 ft was attached to the bottom of the 
package with GE 7031 varnish in order to improve the temperature control, and hence long term stability required for 
this work. The temperature of the vaporizer was controlled with a LakeShore DR-91C Temperature Controller and a 
calibrated CGR-1-2000 carbon-glass sensor, while the sample temperature was controlled with a LakeShore DR-93C 
Temperature Controller and a calibrated CGR-1-2000 carbon glass sensor. In order to achieve long term temperature 
stability, the heater on the vaporizer was set to attain temperature slightly below the desired measurement point; the 
heater on the package was used to fine-tune and stabilize the temperature. Thus, temperature stability of ~ lmK 
could be obtained and held for hours, enabling extremely detailed study of Z s (H)\ Ut T- 

The microwave transmission amplitude S21 of the fully assembled system was measured using an HP 8510C Auto- 
matic Network Analyzer (ANA). All data acquisition and computer control was carried out on an HP 745i workstation 
running HP VEE. 

In order to obtain data at these low field scales it becomes essential to go beyond simply characterizing the resonance 
by measuring the resonance frequency fo and the 3 dB bandwidth A/3 dB and hence the quality factor Q = /0/A/3 dB- 
Thus we developed a least squares analysis to fit our resonances to a Lorentzian form. The results of these fits were 
found to agree very well with the values obtained from the standard 3 dB method, but provided significantly enhanced 
sensitivity to small changes. 

The changes in surface impedance associated with the application of magnetic field are extracted from the measured 
quantities as follows: 

AZ a (H,T)=Z s (H,T)-Z s (H = 0,T) 

= T [Af MB {H, T) - Af 3dB (H = 0,T) + i2 (/ (0, T) - f (H, T))] , 

where T is a geometric factor determined from the sample dimensions. 



2 



III. THEORY 



A. Microwave vortex dynamics 



Viscoelastic models of microwave vortex dynamics, such as those developed early on by Gittleman and Rosenblum 
H and later by Coffey and Clem || consider the effects of the Lorentz force which is exerted on the vortices by the 
microwave currents as the main source of an increase in surface impedance in the presence of magnetic fields. The 
applicability of these models has been validated in earlier experiments on a variety of crystals and films ||-|^] . 

We have shown in Ref. [Q] that the field induced changes in surface impedance of a superconducting thin film can 
be well described using a simple model of vortex dynamics and parameters which are consistent with measurements 
on YBCO single crystals 0, with 

AZ S (H,T) = (T)] B(Hl (2) 

At high fields B(H) —> /j,qH, but this is clearly not the case at low fields or in the presence of bulk pinning where 
the field profile in the sample can be non-uniform. It should also be noted that the increases in surface impedance 
due to penetrated flux must be insensitive to the orientation of the penetrated vortices. Thus, it is insufficient 
to simply replace B{H) with an average flux density B{H) = <&/S = (1/S*) / B(H)dS, where S is the surface of 



the sample perpendicular to the applied field. The desired quantity is the average absolute flux density, \B(H)\ = 
(1/5) J \B(H)\ dS. Thus, AZ S (H,T) can be thought of as counting the absolute number of vortices in the sample 

s 

since it is directly proportional to \B(H)\. Hence for non-uniform flux profiles, the above equation is modified to: 



AZ S {H,T) = \B(H)\, (3) 

fi d[K s (T) + iLur){T)\ 

B. Critical State Flux Penetration Models 

The penetration of magnetic flux into a superconductor depends crucially on two elements, geometry and pinning. 
These are both incorporated in what is usually called the critical state model, originally due to Bean pC[| . Bulk 
pinning is described in terms of the critical current J c , while the field profiles depend heavily on the sample geometry. 
We will compare the results for the more traditional slab geometry in a longitudinal field to those for a thin strip in 
a perpendicular field 

In the following we assume that H c \ is negligible (i.e. that flux penetrates for all fields), and that J C (B) = J c (0) 
(sometimes referred to as the Bean-London model). Although hysteresis loops are usually described in terms of the 
magnetization M(H), the average absolute induced flux density \B(H)\ is the relevant parameter to use in our case, 
as described above. 



1. Bulk Geometries (Slab) 

For a slab of infinite extent of width W, the flux profiles can easily be written down and one can immedi- 

W/2 

ately compute \B{H)\ = (1/W) J \B(x)\dx. When all vortices point in the same direction, \B{H)\ = $/W = 

-W/2 

W/2 

(1/W) J B(x)dx, where $ is the induced flux per unit length in the sample. As the applied field is increased flux 

-W/2 

penetrates gradually from the edges of the slab, until a characteristic field i? s ,buik = J c W/2 is reached and the two 
flux fronts meet at the center of the slab. For increasing field applied to a slab initially in the virgin state 

JBJHY\ = ^ 7 &-^ iorH<H M (4) 

I H — 2^ s > bu lk' f° r H > -Hs, bulk • 
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When the field is reduced after reaching a maximum field H maK > i? S) bulk) the field and current profiles can be 
constructed using a simple construction. 

fl( x , H, if s ,bulk) H max ) = f(x, H max , i? Si buik) — f(x, H max — H, 2ff Sj b u ik), (5) 

where / can be any of the following quantities J, B or M. /j is that same quantity as the field is reduced from H maK . 
Thus, as the applied field is reduced, the flux density at the edges is gradually reduced. At H = the standard 
sandpile flux profile can be found at the center of the sample. As the field is now increased in the other direction, 
negative vortices enter at the edges but the positive vortices at the center remain until \H\ > -ff s ,buik- The curvature 
in the response arises from the insensitivity to direction of flux, and would not be observed in quantities such as 
M(H) or &{H). The calculated hysteresis in \B(H)\ for a slab is shown in Fig. || 



2. Thin Film in a Perpendicular Magnetic Field 

A number of theorists Jl2| fL5| have recently studied the problem of thin strips or disks in a perpendicular magnetic 
field, using conformal mapping techniques introduced early on by Norris |l6) . Analytic solutions for the current and 
flux density distributions in the strip are thus now becoming available for this problem, which could heretofore only 
be considered numerically. 

In what follows we will briefly summarize the central results of Ref. [^2) recast in a form which better suits our 
purposes. Again, as in the case of the traditional Bean model for bulk geometries we ignore H c \ and assume that 
J C (B) = J c (0). 

The geometry is that of an infinitely long strip of thickness d and width 2a aligned such that the sample lies in the 
xy plane with its width along the y axis. In the Meissner state the current distribution must be such that B z (y) = 
for \y\ < a. 

If flux can now penetrate to a distance b < a, the current distribution must change, b can be determined by 
imposing the condition that J(b) is finite, as well as the boundary condition on B {B(x) — (IqH far away from the 
sample. This leads to b = a/ cosh (H/ -ff s , strip), where -ff s , strip is now a characteristic field related to the critical current 
and the geometry, iJ S! strip = Jcd/n. It should be noted that unlike at i? s ,siab in the case of a slab, the two flux fronts 
do not meet at H = i? S! strip- In fact, b is always non-zero for finite fields, although at some point b will be smaller 
than the effective size of a single vortex A and the validity of these expressions must be questioned. The resulting 
sheet current distribution can then be written as 

K x (y)^ I J x {z,y) = { * (6) 

Jcd-fo, b < \y\ < a. 

'0, \y\<b 
BM= ^ogctanh- 1 gggr , b < \y\ < a (7) 
M^tanh- 1 ^^, o<|y| 

It should be noted that B z (y) has a logarithmic singularity at the edges in this expression. 
Thus, we can use Eq. |?] to evaluate 

a 

\BjH)\ = (l/2a) f B(y)dy 

(8) 

TJ 

= Mo# S) stri P In cosh — 

**-s, strip 

Note that this has some striking similarities with the results for the slab geometry: (i) a characteristic field scale 
(-ff^buik or -ff S! stri P ) related to the critical current density, (ii) \B(H)\ cx H 2 at low fields H <C H s , and (iii) \B(H)\ cx H 
at high fields H > H s . 

Assuming that Eq. || applies just as for the slab, the hysteretic response can be constructed with -ff S! strip taking the 
place of -ff^buik, and is presented in Fig. |3|. 
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IV. Z S {H): EXPERIMENTAL RESULTS 



All of the results presented here were obtained with the sample at fixed temperature, and initially in the virgin 
state. This allows us to separate out the initial penetration of flux from the overall response and hysteresis which is 
observed at higher fields. The sample was warmed up well above T c and cooled down in zero field between runs in 
order to ensure that the sample was really in the virgin state. 

The data for the virgin state flux penetration from < 1 Oe to 6 T fall naturally into 3 field regimes: 

1. changes occur initially on very small field scales which we call Hi(4.2K) ~ lOOOe. 

2. This is followed by the majority of the flux penetration at field scales of around H S (A.2K) ~ 5000 Oe character- 
ized by a superlinear dependence of the impedance on H, and indicating an inhomogeneous flux distribution. 

3. At very high fields H H s the dependence of Z S (H) becomes linear in H indicating uniform flux penetration. 



A. Initial Penetration of Flux 

The detailed behavior up to 400 Oe at A.2K is presented in Fig. [|, where it can clearly be seen to be characterized 
by a field scale Hi(4.2K) ~ lOOOe, and that above and below Hi the impedance appears to be approximately linear. 
The actual dependence is well described by the functional form: 

\B(H)\ = a taah(H/Hi) + QH (9) 

and is shown in the figure as a solid line. The value of Hi used in the fit is indicated by an arrow. 

It is easy to extract Hi from the data by a linear extrapolation of the behavior above and below the bend and is 
presented in Fig. ^| Hi is chosen to be the intersection point of these two lines. This is nearly equivalent to the field 
scale Hi that enters in Eq. g, but does not presuppose any specific functional dependence of Z S (H) at low fields. 

Is this low field scale simply a manifestation of H c \l H c \ is typically observed as a sharp break from practically zero 
in microwave measurements, due to the logarithmic singularity in B(H) near H c \. This is not what is seen in Fig. |J. 
However, the magnetic field i? c i,strip a t which flux initially penetrates into a thin strip in a perpendicular magnetic 
field will in general be lower than the bulk H c \ of the material, due to demagnetization effects. Using standard 
ellipsoidal approximations to a thin strip yields an estimate of -Hci, strip = H c \{dj2a). Recent analytical calculations 
find that the effect is actually smaller than that predicted by the ellipsoidal approximation by a factor of y d/2a, 
so that H c i strip = H c iy / d/2a. For a 0.5 /zm x 100 /im strip, typical of those used for these experiments, one finds 

that i?cl : strip ~ flcl/10. 

Bulk H c i for single crystalline YBCO was measured by Wu et al. where they found that H c i ,||(0) = 850 ±40 Oe, 
leaving us with an upper bound on the penetration field -ff c i, strip ~ 85 ± 4 Oe for our samples. This magnitude is 
comparable to the observed low temperature magnitude of Hi. 

If Hi is simply a manifestation of H c \, the linear field dependence below Hi is somewhat unconventional. For 
a conventional GL superconductor, the effects of pairbreaking are expected to have a quadratic field dependence 
below H c i j^),^. This quadratic dependence has also been observed experimentally in radio frequency penetration 
depth measurements on polycrystalline samples an d single crystals of YBCO for fields applied parallel the the 
ab plane. It should also be mentioned that Yip and Sauls pl[ have shown that for a d -wave superconductor, the 
penetration depth (and hence also the surface resistance) should rise linearly in the Meissner state. 



1. H c i and Weak Links 

A quadratic field dependence is also expected from the weakly coupled-grain model j2^j2j|, which has been used 
to explain the functional dependence of Z s on low microwave fields. Certain theories for granular superconductors do 
predict an initial linear field dependence (24|,|25|] by considering a Bean-like model applied to the inter-granular weak 
links. Weak links are proposed as an explanation for a similar initial linear rise followed by a saturation reported in 
the literature pfifl . We looked for, but did not find any direct correlation between film granularity and the magnitude 
Hi. Films of very different granularity (as evidenced by the rate at which the films etched) had remarkably similar 
magnitudes of Hi. Thus, the linear dependence on field cannot conclusively be attributed to the presence of weak 
links in the samples. 
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B. Scaling of the data for moderate fields 



The field induced change in surface impedance, AZ s (H,T,ui) = Z s (H,T,uj) — Z S (H — > 0, T, 10) at moderate fields 
(Hi < H < 4 T) scales, giving a field scale H s and an impedance scale which can be determined from the slope at 
high fields. 

The scaled data over this broad field range are well described by the following expression 



\B(H)\ = C s \ncosh{H/H s ) (10) 

where H s is the field scale well above which the flux has largely penetrated throughout the sample and B(H) « {IqH . 
The scaling function used above is exactly that obtained from the critical state model of rectangular strips presented 



in Sec. IIIB 2. This would imply that J c must depend only weakly on H, and only one field scale (related to J C (T)) 
determines the response. Thus, we should now be able to estimate a critical current density J c from the experimentally 
observed field scale and the film thickness using H s = J c d/%. 

A similar scaling was also observed in data taken from a ring resonator. Calculations of \B(H)\ for the case a 
superconducting disk using the results of Ref. |D| cannot be distinguished from the lncosh(_ff/iJ s s t r i p ) behavior of 
the strip within experimental uncertainty, although the relation of the scaling field -ffs.disk to J c is different. We 
have presented the temperature dependence of J C {T) in Fig. ^ for the ring resonator, obtained using Eq.|lO|. (This 
overestimates the true J c values by a factor of up to 7r/2, due to the ring geometry although the magnitudes are 
consistent). 

When combined with the J c values estimated from the nonlinear response at high T on the same resonator 
from Ref. || it would appear that the two estimates are consistent, and can well be described overall by a 
(1 — t 2 ) 2 dependence reminiscent of a mean field thermodynamic critical field, H c . Also the magnitudes obtained, 
J c (77 K) ss 2 x 10 6 A/cm 2 , comparable with independent DC measurements on films of the same batch. 

To summarize, flux penetration into the virgin state for moderate fields is very well characterized quantitatively by 
the critical state model for perpendicular geometry. 



C. High field response for H ^> H s 



For very high fields much greater than the self field H s the impedance becomes linear in the applied field. In this 
regime the slopes can be used to extract vortex parameters such as the pinning force constants and the viscosity as 
we have done in Ref. 101. 



V. HYSTERESIS 



So far we have only examined what happens as the externally applied magnetic field is increased adiabatically. The 
response as the field is slowly reduced and also reversed will be discussed in this section. 



A. Moderate Field Hysteresis 



The typical response when the applied field is ramped up to a field significantly higher than H s and then reduced 
is presented in Fig. ||. As can immediately be seen this response is qualitatively different from that predicted by the 
critical state model alone (presented in Fig. |J). The loop has the opposite handedness, and is much smaller than H s 
in width. While the reduced width of the loop may be consistent with the presence of a geometrical or other barrier 
at the strip edge, the handedness however cannot be so easily explained. 

There are a number of possible ways to explain this behavior: 

• The relationship between Z s and \B(H)\ does not apply, 

• The construction (Eq. |^) used to determine B± (H) does not apply, or 

• A fraction of the flux, principally at the edges, is sufficiently weakly pinned to be able to leave the sample 
quickly as the applied field is reduced. 
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B. Hysteresis at Very Low Field 



In order to shed some light on the issue, we have also carried out measurements of the hysteretic response for fields 
~ Hi. Typical response when the applied field is ramped up to a field i? m ax comparable to Hi and then down to 
-fl max and finally increased back to H max is presented in Fig. 0. Although only AX S (H) is presented in the figure, 
AR S (H) has qualitatively the same behavior. At sufficiently high T, response qualitatively similar to that observed 
in the previous section is recovered even for fields \H\ < 400 Oe. 

To put these effects into perspective Fig. [ll] presents the data of Fig. ||, along with low field data obtained at the 
same temperature. Note that the magnitude of these effects is very small, and requires extreme care to be observed. 

Clearly the clockwise nature of the hysteresis initially observed at high fields is preserved even at these low field 
scales, and the drop in Z s as the applied field is reduced appears to be rather sudden, taking place within < 10 Oe. 

Below we summarize the essential features of the hysteretic response: 

• Sudden drop in both R s and X s as the applied field is reduced, below even the virgin state response 

• After -ffmax is reached, the response is clearly symmetric in H about H = 0. 

• Z s at H = does not return to its initial value in the normal state. 

• a cusp-like feature is observed as H is increased in the opposite direction after reaching H = 0. 

A cusp-like feature is expected in magnetization as the field crosses H — in the presence of a geometrical barrier 
. One must also keep in mind that the time taken to manually switch the current leads (to reverse the field) was 
finite in these experiments so that any time dependence (such as flux creep) may also be a possible source for the 
observed cusp. 

C. Comparison to DC Magnetization 

Hysteresis loops from DC magnetization measurements on thin films and bulk polycrystalline samples in high fields 
(H >^ I T) generally appear quite different from those we report here, and are often more consistent with critical 
state models. M(H) loops are generally counter clockwise and estimates of J C (H) are routinely extracted from the 
directly from the loop width. Measurements on single crystalline samples also have the "correct" handedness although 
the fishtail effect is observed in the vast majority of samples. 

This is in contrast to the results presented here where the hysteresis loops have the opposite handedness are much 
narrower than expected from the critical state model. In fact, rather than being related to H s , the hysteresis width 
appears to be governed by a field scale on the order of Hi. 

VI. CONCLUSION 

The linear electrodynamic response of vortices in thin films in the perpendicular geometry has been investigated 
within the framework of recent theoretical work. While flux entry into a sample initially in the virgin state appears 
to be well described by the critical state model, flux exit when the field is reversed does not follow the model. The 
expected hysteresis from the critical state model ( of magn itude « 2H S ), is not observed. In fact, the observed 
behavior would imply that the absolute flux density \B(H)\ drops suddenly as the applied magnetic field is reduced. 
The presence of geometrical barriers at the strip edges can help to reduce the level of hysteresis expected, but will 
not lead to the observed response. 
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Resonator 


Geometry 


Line Width 


Resonator Length 


Film Thickness 


Substrate Thickness 


N4 


Straight 


100 /im 


1.32 cm 


0.38 /im 


250 fim 


N3 


Ring- 


100 ^m 


2.00 cm 


0.375 fim 


500 Atm 



TABLE I. The two resonators used for this work 



FIG. 1. Typical virgin state response: Raw A/ data for increasing applied field for a sample initially in the virgin state for 
resonator N3 at 3.7 GHz and 5 K. The dashed line is the expression of Eq. |l^. The arrows indicate the two field scales, Hi and 
H s . 



FIG. 2. The calculated hysteresis in \B(H)\ for a thick slab in a field applied parallel to its surface. Both axes are in units 

Of -ffs.bulk, with ff maJ = lO-ff^bulk. 



FIG. 3. Hysteretic response in |_B(i/)| for a thin strip in a perpendicular field calculated numerically from the field profiles 
discussed in the text. Both axes are in units of -Hs, strip, and // max = 10// a , strip- 



FIG. 4. Typical virgin state response: Raw A/ data for increasing applied field for a sample initially in the virgin state for 
resonator N4 at 3.7 GHz and 5 K. The dashed line is the expression of Eq. ^. The arrow indicates the field scale described in 
the text. 
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FIG. 5. Temperature dependence of the initial penetration field scale for resonator N4. The straight dashed line is a least 
squares fit to the data. The dotted line presents a BCS temperature dependence for H c \. 



FIG. 6. Raw (left) and scaled (right) plot of the field induced changes is surface resistance. Data presented is for resonator 
N3 at 3.7 GHz and T = (5.0, 15.0, 20.0, 25.0, 31.0, 40.0 and 50.0 K). The slope increases with increasing T temperatures. 
The solid line represents the expected behavior for a strip. A dashed line representing the expected behavior for a disk is also 
presented, but is indistinguishable from the solid line. 



FIG. 7. Field scale H s for resonator N3 at 3.7 GHz. The dashed and dotted lines are (1 — t 2 ) and (1 — t 2 ) 2 dependences 
respectively 



FIG. 8. Hysteretic response of R S (H) for ff ma x 3> H s . Results are presented for resonator N4 at 3.7 GHz and 10 K. 



FIG. 9. Hysteretic response of X a (H) for // max ~ Hi. Results are presented for resonator N4 at 3.7 GHz and 35 K. 



FIG. 10. Low field hysteresis in R S (H) at high temperatures. Results are presented for resonator N4 at 3.7 GHz and 84 K. 



FIG. 11. Magnified view of the hysteresis data for R s presented in Fig. g for moderate fields and low fields (inset). 
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